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On the nature of the residual meson-meson interaction from simulations 
with a QED2+1 model * 
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A potential between mesons is extracted from 4-point functions within lattice gauge theory taking 2+1 dimen- 
^ sional QED as an example. This theory possesses confinement and dynamical fermions. The resulting meson- 
meson potential has a short-ranged hard repulsive core and the expected dipole-dipole forces lead to attraction 
I at intermediate distances. Sea quarks lead to a softer form of the total potential. 
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1. Introduction. The extraction of an effec- 
tive interaction, or potential, between two com- 
posite hadrons from a quantum field theory is of 
fundamental importance. Here we have in mind 
the quark-gluon degrees of freedom at the sub- 
nuclear level. Those are thought to give rise to 
the nucleon-nucleon forces usually described by 
meson theory leading to the construction of the 
Bonn and Paris potentials [1]. 

In the last years lattice calculations with static 
quarks have demonstrated that the potential bet- 
ween two three-quark clusters is attractive [2]. 
A hard repulsive core of the potential, as sug- 
gested by experiments and their interpretation, 
could not be observed in the region where the 
two nucleons have relative distance close to zero. 
Beyond the static theory, employing dynamical 
quark propagators, correct antisymmetrisation 
and quark exchange become possible [3,4]. In this 
framework the mechanism which leads to a repul- 
sive core will be investigated here analytically and 
numerically within lattice QED2+1. 

2. Four-Pomt Correlation Functions. We 
take the one-meson field as a product of staggered 
Grassmann fields x ^nd x with external fiavours 
u and d 

'i>x{t) = Xd{xt)xu{xt) . (1) 
The meson-meson fields with relative distance r = 



y — X are then constructed by 



(2) 



The information about the dynamics of the 
meson-meson system is contained in the 4-point 
time correlation matrix 

C^,{tM = (<I>tW<I>f'(^o))-(<I>tW)(<I>f'(^o)).(3) 

Working out the contractions between the Grass- 
mann fields the following diagrammatic represen- 
tation is obtained 



^,(4A) _^ (j(AB) _ (j(AC) _ (j(AD) (-4^ 



+ X - M -M ^ w 

Each of the four contributions to the correlator 
comprises the exchange of gluons and sea quarks. 
For diagrams C(4^) and C(4^) those take place 
between the mesons, whereas diagrams and 

C(AD) 

correspond also to interaction mediated by 
the exchange of valence quarks. Denoting the 
contracted Grassmann fields by the quark propa- 
gators Gn, n = 1 . . .4, we have for example 



(j(AA) 
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(6) 
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y ^ X ^ X ^ y I 

where ~ stands for the sums with normalisation 
factors that carry over from (2). 

In order to define an effective meson-meson in- 
teraction it is crucial that the noninteracting com- 
ponents in are isolated [4]. To the corre- 
sponding expression, = C''-'*^-' -|- &^^\ only 
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the direct channels contribute. This respects the 
boson symmetry on the meson-field level and con- 
stitutes the free meson-meson correlator with the 
relative interaction switched oS. 

3. Effective Interaction from Time Evolution. 
The deviation of C'-'*-' from contains the 

residual effective meson-meson interaction. From 
analogy with quantum mechanical perturbation 
theory 



C{t-to) = e-*^(*-*°) 



(7) 



we define an effective meson-meson interaction as 



■Hi 



dC 

dt 



,withC = C'(4)""C7(4)C'(4)"". (8) 



This definition is meant as a relation between ma- 
trices and holds independently of the basis. It 
can be tested both in momentum space [5] and 
in coordinate space. The above equations build a 
bridge between the quantum field theoretical cor- 
relation functions on the lattice and an effective 
quantum mechanical Hamiltonian. 

Jf. Hard Core from Adiahatic Approximation. 
An analysis of the coordinate-space matrix el- 
ements of the correlators C'-'*-' and reveals 
an interesting perspective on the repulsive core. 
From C'-'*-' expressed in terms of the fermion prop- 
agators it is easy to see analytically that 



C%\{t,to) = 0, if r = Oorr' = 0, 



(9) 



In our simulation we choose the relative distance 
between the mesons to be the same at the initial 
and final times of the propagation, r = r' . In the 
spirit of the Born-Oppenheimer approximation it 
is assumed that the interaction proceeds faster 
than the motion between the mesons. This seems 
to be justified because the quark mass mp = 0.1 
is relatively large. In this sense we replace the 
sum over all eigenstates in the spectral represen- 
tation of the correlation function by an average 
term with an effective energy W{r) 



r|n)|2e-^''(*-*°) 



4?(Mo) = E 

n 

~ c(r)e-^('^)(*-*°) 



(10) 



To the extent that the energy W{r) of the meson- 
meson system at fixed relative distance r can be 



extracted from the large-t behaviour of the diago- 
nal elements we may conclude from (9)-(10) that 
W{r = 0) = -|-cxD, provided the strength factor 
c(r ) is nonzero at r = 0. In this case the effec- 
tive interaction possesses a hard repulsive core. 
It is due to the anticommuting nature of the con- 
stituent fermion fields leading to Pauli repulsion. 

5. Results. The simulation was performed 
on an 1,2 X T = 242 X 32 periodic lattice with 
the compact ?7(1) Wilson action at /3 = 1.5 and 
staggered fermions with fiavour number Np = 
and Np = 2, both with dynamical quark mass 
mp = 0.1. We used 64 independent gauge field 
configurations generated with the molecular dy- 
namics algorithm [6]. Fermion propagators were 
computed with 32 and 16 random sources [7] in 
the quenched case and in full QED2+1, respec- 
tively. 

The potentials were extracted from cosh-fits to 
the correlator ratios C^^^^ for t = 15 ... 19 in 

periodic time. The first definition (7)-(8) repre- 
sents a matrix and leads to a nonlocal Hamil- 
ton operator 7ii(f,f') where r,f' is the sepa- 
ration between the mesons. Since computation 
of the entire 4-point correlation matrix is very 
time-consuming we only obtained the diagonal el- 
ements. This also allows a direct comparison with 
>V/(r) = W(r) — 2m from the second definition 
(10), where 2m is the mass of the noninteracting 
two-meson system. 

We begin the description of the results of our 
simulation with the quenched theory, Np = 0. 
In fig. 1 we display data of the correlation func- 
tion Cff for the direct channels (A-|-B), the quark- 
exchange channels (C-|-D) and for their combina- 
tion (A-|-B)-(C-|-D). One observes a zig-zag be- 
haviour for even and odd distances r. In the 
staggered scheme the meson correlator exhibits 
an oscillating term cx (— l)*e~'" * with m' almost 
degenerate with the ground state, m' m. On a 
euclidean lattice there is a corresponding term in 
the space directions. This feature enters into the 
meson-meson correlator. For on-axis points r the 
large-t behaviour of the correlator is 



« c(+)(r)e-^i C-^)* + (_)'-c(-)(r)e-^r'«*. 

If one is interested in the meson-meson inter- 
action between ground states only it is neces- 



Figure 1. Raw data of correlator ratios C (sym- 
bols) and coefficients c^^) from an even-odd anal- 
ysis as a function of relative distance r at t = 17. 
Comparison of the (a) direct, (b) quark-exchange 
and the (c) combined channels. 



sary to separate c^^). Independent 3-parameter 
fits to even and odd sets of r were made using 
/(r) = aie~"^'' + a^. Remarkably, the pairs of 
interpolated curves c^^) in fig. 1 exhibit very sim- 
ilar shapes (aside from different clustering). This 
indicates that the m-m and m'-m' interactions 
are similar, thus justifying the use of the stag- 
gered scheme. 

By repeating the above procedure for fixed t 
and subsequent cosh-fits in t-directon we recon- 
structed the correlator ratios for the ground state 
meson-meson system. The results for Ti^'j^^r) 

are presented in fig. 2a and those for Wj+V) 
in fig. 2b. The (— ) data are lost in noise. Both 
direct diagrams (A-|-B) correspond to a repulsive 
potential barier. The quark-exchange diagrams 
(C-|-D) yield similar behaviour (not shown) ex- 
cept their contributions enter with a minus sign 
into the total correlations (4), (5). This gives an 
attractive total potential as displayed in fig. 2. 

We now consider correlations for all available 
distances r on the axes and diagonals. In fig. 3 the 
corresponding potentials of the direct diagrams 
(A-l-B) and all four diagrams (A-|-B)-(C-|-D) are 
shown. The direct channels turn out to be repul- 
sive as in fig. 2. For both total potentials %i{r) 
and yVi{r) there is now some evidence for a repul- 
sive core, as discussed earlier, followed by attrac- 
tion at medium distances. The shapes emerge as 
a subtle interplay between all graphs (A-|-B) and 
(C+D). 
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Figure 2. Meson-meson potentials for the di- 
rect channels (A-|-B) and all four graphs (A-|-B)- 
(C-|-D) extracted from interpolated correlators. 
The result (a) 7ii(r) from the effective time evolu- 
tion is compared with (b) >V/(r) of the adiabatic 
approximation. In both cases the direct channels 
lead to a repulsive potential barier. The total 
potential comes out attractive. Curves are guide 
to the eye and error bars result from a jack-knife 
analysis. 



The Pauli principle leads us to expect a hard 
repulsive core. Our data show that its range is 
short. This is in contrast to a complementary 
simulation in a truncated momentum basis [4,5]. 
Therefore, the actual dynamical size of one meson 
on the lattice is of interest. Consider the fields 

<l'p{^) = ^ 6p,x-yXd{yt)Up{yt, xt)xu{xt) , 

where the product of link variables U^{yt,xt) 
along the quark-antiquark separation p ensures 
gauge invariance of the extended meson operator. 
The corresponding correlation matrix ii{t,to) 
was computed for p, p' = 0, 1,2. Reduced matri- 
ces of C*-^-* for all relevant irreducible representa- 
tions of the lattice symmetry group 0(2, Z) were 
obtained. It turned out that the matrix elements 
with Pj p' z/z were numerically zero. This indi- 
cates that, for the current model, the meson ac- 
tually stays small (< a) during its propagation. 
With regard to the meson-meson interaction, this 
explains the narrow hard core. 

Simulations for Np = 2 fiavours also have 
been performed. Using only 16 random sources 
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Figure 3. Meson-meson potentials as in fig. 2, 
but from raw correlators. The direct channels are 
repulsive again. Here both total potentials hint 
at a repulsive core at short range. 



Figure 4. Meson-meson potentials as in fig. 2, 
extracted from full QED2+1. The direct channels 
(open symbols without error bars) are small. The 
total potentials resemble those of the pure gauge 
case. 



the data are considerably more noisy. For 
ease of comparison, we rescale in fig. 4 Tij 
and W/ by a factor equal to the mass ratio 
2m(NF = 0)/2m(NF = 2) = 1.031(9)/0.63(1) = 
1.65. For the total potentials there seem no qual- 
itative changes compared to the pure gluonic case, 
except that sea quarks apparently lead to a softer 
behaviour. The shapes of W/ again resemble Tij . 
The potentials from the direct channels are small. 

6. Conclusion. We have proposed two pos- 
sibilities to extract an effective interaction bet- 
ween two composite particles. We found for 
QED2+1 that gluon exchange and valence-quark 
exchange are equally important. The definite 
shape of the potential might depend on sea quark 
effects. In the model used mesons feel a short- 
ranged hard-core potential due the Pauli princi- 
ple. The medium-range residual forces between 
colour-neutral quark-clusters are attractive. It 
should be noted that a computation of meson- 
meson scattering phase shifts based on Liischer's 
proposal [8] indicates the same behaviour of the 
potential [9] as do a momentum-space study [4,5] 
and a different coordinate-space approach [3]. 
The correct antisymmetrisation of the meson- 
meson correlator was important for the emer- 
gence of a hard core. In contrast to U(l), the 
nonabelian colour structure of SU (3) will not nec- 
essarily lead to a vanishing correlator for nucleon- 
nucleon distance zero. Preliminary studies seem 



to indicate attraction. Calculations with light va- 
lence quarks might help to understand the hard- 
core problem of the nucleon-nucleon interaction. 
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